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Resumo 

Quando cercado por uma região de emissão trans¬ 
parente, espera-se que os buracos negros revelem uma 
sombra escura causada pelo efeito da gravidade na 
curvatura da luz e pela captura de fótons no horizonte 
de eventos. Para ilustrar e estudar este fenômeno, 
montamos o Event Horizon Telescope (EHT), um ar¬ 
ranjo interferométrico de ampla linha de base obser¬ 
vado no comprimento de onda de 1,3 mm. Isso nos 
permite reconstruir imagens em escala do horizonte 
de eventos do candidato a buraco negro supermassivo 
no centro da galáxia elíptica gigante M87. A fonte 
de rádio compacta central foi resolvida como um anel 
de emissão brilhante assimétrica com um diâmetro de 
42±3/ras, circular e englobando uma depressão central 
com brilho em uma razão de fluxo > 10 : 1. O anel 
de emissão é recuperado usando diferentes esquemas 
de calibração e imagem, com seu diâmetro e espessura 
permanecendo estáveis ao longo de quatro observações 
realizadas em diferentes dias. No geral, a imagem 
observada é consistente com as expectativas para a 
sombra de um buraco negro de Kerr, como previsto 
pela relatividade geral (RG). A assimetria no brilho 
do anel pode ser explicada em termos de irradiação re- 
lativística da emissão de um plasma girando próximo 
da velocidade da luz ao redor de um buraco negro. 
Comparamos nossas imagens a uma extensa biblioteca 
de simulações magneto-hidrodinâmicas de buracos ne¬ 
gros usando a RG e derivamos uma massa central de 
M = 6,5 ± 0,7 X 10®M©. Nossas observações em 
ondas de rádio fornecem portanto evidências podero¬ 
sas sobre a presença de buracos negros supermassivos 
em centros de galáxias e como motores centrais de 
núcleos ativos de galáxias. Eles também apresentam 
uma nova ferramenta para explorar a gravidade em 
seu limite mais extremo e em uma escala de massa 
até agora inacessível. 


*Este artigo foi traduzido por Sthiven Raphael Melo 
Correia entre lO/abril/2019 e 31/maio/2019 com o auxilio 
do pós-doutorando do INCT do e-Universo (LIneA/ON) 
Adriano Pieres. A intenção é prover maior acessibili¬ 
dade do artigo aos alunos que apresentam dificuldades na 
leitura no texto original (https://iopscience.iop.org/ 
article/10.3847/2041-8213/ab0ec7/pdf). 


1 Introdução 


Os buracos negros são uma previsão funda¬ 
mental da teoria da relatividade geral (RG; Eins- 
tein 1915). Uma característica principal dos bu¬ 
racos negros é seu horizonte de eventos, um limite 
causal unidirecional no espaço-tempo do qual 
nem mesmo a luz pode escapar (Schwarzschild, 
1916). A produção de buracos negros é genérica 
em RG (Penrose 1965), e mais de um século de¬ 
pois de Schwarzschild, eles permanecem no cen¬ 
tro das questões fundamentais na unificação da 
RG com a física quântica (Hawking, 1976; Gid- 
dings, 2017). 

Buracos negros são comuns em astrofísica e 
são encontrados em uma ampla gama de mas¬ 
sas. Evidência de buracos negros de massa este¬ 
lar vem de raios-X (Webster & Murdin 1972; Re- 
millard &: McClintock 2006) e medições de ondas 
gravitacionais (Abbott et al. 2016). Considera- 
se que buracos negros supermassivos, com mas¬ 
sas de milhões a dezenas de bilhões de massas 
solares, existem nos centros de quase todas as 
galáxias (Lynden-Bell, 1969; Kormendy e Ri- 
chstone, 1995; Miyoshi et ah, 1995), inclusive 
no centro da nossa Galáxia (Eckart & Genzel 
1997; Ghez et al. 1998; Gravity Collaboration 
et al. 2018a) e no núcleo da galáxia elíptica M87 
(Gebhardt et al. 2011; Walsh et al. 2013). 

Os núcleos galácticos ativos (AGNs) são 
regiões brilhantes centrais que podem ofuscar 
toda a população estelar de sua galáxia hospe¬ 
deira. Alguns desses objetos, os quasares, são as 
fontes mais luminosas do universo (Schmidt 1963; 
Sanders et al. 1989) e se considera que são ali¬ 
mentados por buracos negros supermassivos que 
acretam matéria a taxas muito altas através de 
um disco geometricamente fino e opticamente es¬ 
pesso (Shakura & Sunyaev 1973; Sun & Mal- 
kan 1989). Em contraste, a maioria dos AGNs 
no universo local, incluindo o centro galáctico e 
M87, estão associados a buracos negros super¬ 
massivos alimentados por quentes e débeis fluxos 
de acreção, com taxas de acreção muito menores 
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(Ichimaru 1977; Narayan e Yi 1995; Blandford e 
Begelman 1999; Yuan & Narayan 2014). 

Em muitos AGNs, os jatos de plasmas rela- 
tivísticos colimados (Bridle & Perley 1984; Zen- 
sus 1997) lançados pelo buraco negro central con¬ 
tribuem para a emissão observada. Esses ja¬ 
tos podem ser alimentados tanto por campos 
magnéticos que segregam o horizonte de eventos, 
extraindo a energia rotacional do buraco negro 
(Blandford & Znajek, 1977), quanto pelo fluxo 
de acreção (Blandford e Payne, 1982). A emissão 
no horizonte próximo dos núcleos galácticos ati¬ 
vos de baixa luminosidade (LLAGNs; Ho, 1999) é 
produzida pela radiação síncrotron, que atinge o 
pico do rádio através do infravermelho distante. 
Esta emissão pode ser produzida pelo fluxo de 
acreção (Narayan et al. 1995), pelo jato (Ealcke 
et al. 1993), ou por ambos (Yuan et al. 2002). 

Quando visto do infinito, um buraco negro não 
rotativo de Schwarzschild (1916) tem um raio de 
captura de fótons Rc = onde rg = GM/c^ 

é a escala de comprimento característica de um 
buraco negro. O raio de captura de fótons é 
maior que o raio de Schwarzschild Rs que marca 
o horizonte de eventos de um buraco negro não 
rotatório, Rs = 2rg. Eótons se aproximando 
do buraco negro com um parâmetro de impacto 
b < Rc são capturados e mergulham no buraco 
negro (Hilbert 1917); fótons com b > Rc escapam 
ao infinito; fótons com b = Rc são capturados 
em uma órbita circular instável e produzem uma 
lente comumente chamado de “anel de fótons”. 
Na métrica Kerr (1963), que descreve buracos ne¬ 
gros com momento angular de giro, Rc muda com 
a orientação do raio em relação ao vetor momento 
angular, e a seção transversal do buraco negro 
não é necessariamente circular (Bardeen 1973). 
Essa mudança é pequena (< 4%), mas potencial¬ 
mente detectável (Takahashi 2004; Johannsen &: 
Psaltis 2010). 

As simulações de Luminet (1979) mostraram 
que para um buraco negro imerso em um disco de 
acreção geometricamente fino e opticamente es¬ 
pesso, o raio de captura de fótons pareceria a um 
observador distante como um fino anel de emissão 
dentro de uma imagem do disco de acreção seme¬ 
lhante a uma lente. Para buracos negros acreta- 
dores imersos em um disco opticamente fina e 
geometricamente espessa, como nos LLAGNs, a 
combinação de um horizonte de eventos com a 
curvatura da luz leva ao aparecimento de uma 
“sombra” escura acompanhado de um anel de 
emissão brilhante que deve ser detectável através 
de experimentos de interferometria de linha de 
base ampla (Ealcke et al. 2000a) ou VLBI, na 
sigla em inglês. Sua forma pode aparecer como 


um “crescente” por causa da rotação rápida e ir¬ 
radiação relativística (Ealcke et al. 2000b; Brom- 
ley et al. 2001; Noble et al. 2007; Broderick & 
Loeb 2009; Kamruddin Sz Dexter 2013; Lu et al. 
2014). 

O diâmetro projetado observado do anel de 
emissão, que contém radiação principalmente do 
anel de fótons gravitacionalmente em forma de 
lente, é proporcional ao Rc e, portanto, à massa 
do buraco negro, mas também depende não tri¬ 
vialmente de vários fatores: a resolução de ob¬ 
servação, o vetor rotação do buraco negro e a sua 
inclinação, bem como o tamanho e a estrutura 
da região de emissão. Esses fatores são tipica¬ 
mente da ordem da unidade e podem ser calibra¬ 
dos usando modelos teóricos. 

Simulações modernas de fluxos de acreção 
usando RG e transferência radiativa produzem 
imagens realistas de sombras e crescentes de bu¬ 
racos negros para uma ampla gama de modelos 
de emissão próximos ao horizonte (Broderick & 
Loeb 2006; Moscibrodzka et al. 2009; Dexter et 
al 2012; Dibi et al. 2012; Ghan et al. 2015; 
Moscibrodzka et al. 2016; Porth et al., 2017; 
Chael et al., 2018a; Ryan et ah, 2018; Davelaar et 
al., 2019). Estas imagens podem ser usadas para 
testar as propriedades básicas dos buracos negros 
como previsto em RG (Johannsen & Psaltis 2010; 
Broderick et al. 2014; Psaltis et al. 2015), ou em 
teorias alternativas da gravidade (Grenzebach et 
al. 2014; Younsi et al. 2016; Mizuno et al. 2018). 
Eles também podem ser usados para testar teo¬ 
rias alternativas para buracos negros (Bambi & 
Freese 2009; Vincent et al. 2016; Olivares et al. 
2019). 

VLBI em um comprimento de onda de ob¬ 
servação de 1,3 mm (230 GHz) com linhas de 
base da escala de diâmetro da Terra é necessária 
para resolver as sombras do núcleo de M87 (M87* 
daqui em diante) e do centro Galáctico de Sa¬ 
gitário A* (Sgr A*, Balick & Brown 1974), os 
dois buracos negros supermassivos com os mai¬ 
ores tamanhos aparentes angulares (Johannsen 
et al. 2012). Com comprimentos de onda de 
1,3 mm e mais curtos, linhas de base VLBI com 
diâmetro da Terra alcançam uma resolução an¬ 
gular suficiente para resolver a sombra de ambas 
as fontes, enquanto os espectros de ambas as fon¬ 
tes se tornam opticamente finas, então revelando 
a estrutura da região de emissão mais interna. 
Experiências iniciais (Padin et al. 1990; Krich- 
baum et al. 1998) demonstraram a viabilidade 
de técnicas de VLBI em comprimentos de onda 
de ~ 1,3 mm. Ao longo da década seguinte, um 
programa para melhorar a sensibilidade de 1,3 
mm-VLBI através do desenvolvimento de instru- 



mentação de banda larga levou à detecção de es¬ 
truturas com escala de horizonte de evento em 
Sgr A* e M87* (Doeleman et al. 2008, 2012). 
Com base nessas observações, a colaboração Te¬ 
lescópio de Horizonte de Eventos (EHT) foi es¬ 
tabelecida para montar um arranjo global VLBI 
operando em um comprimento de onda de 1,3 
mm com tal resolução angular, sensibilidade e 
cobertura de linha de bases necessárias para ima- 
gear as sombras em M87* e Sgr A*. 

Neste artigo, apresentamos e discutimos as 
primeiras imagens em escala do horizonte de 
eventos do candidato a buraco negro supermas- 
sivo M87* de uma campanha EHT VLBI reali¬ 
zada em abril de 2017 em um comprimento de 
onda de 1,3 mm. Os artigos anexos fornecem 
uma descrição mais detalhada do instrumento 
(EHT Collaboration et al. 2019a, Documento 
H), redução de dados (EHT Collaboration et al. 
2019b, daqui em diante Documento Hl), ima¬ 
gens da sombra M87 (EHT Collaboration et al. 
2019c), daqui em diante, Documento IV), mo¬ 
delos teóricos (EHT Collaboration et ah, 2019d, 
daqui em diante. Documento V), e a estimativa 
da massa do buraco negro (EHT Collaboration 
et al. 2019e, daqui em diante Documento VI). 

2 O Núcleo de M87 em ondas de rádio 


Em Curtis (1918), Heber Curtis detectou uma 
característica linear em M87, mais tarde cha¬ 
mada de “jato” por Baade &: Minkowski (1954). 
O jato é visto como uma fonte de rádio brilhante, 
Virgo A ou 3C 274 (Bolton et al. 1949; Kassim et 
al. 1993; Owen et al. 2000), que se estende até 65 
kpc com uma idade estimada em cerca de 40 Myr 
e potência ~ 10^^ a 10^^ erg s“^ (de Gasperin et 
al. 2012; Broderick et al. 2015). Também é bem 
estudado na faixa do ótico (Biretta et al. 1999; 
Perlman et al. 2011), Raio-X (Marshall et al. 
2002) e bandas de raios gama (Abramowski et al. 
2012). A extremidade a montante do jato é mar¬ 
cada por uma fonte de rádio compacta (Cohen 
et al. 1969). Tal fonte de rádio compacta está 
presentes em todos os LLAGNs (Wrobel & Hees- 
chen 1984; Nagar et al. 2005) e acredita-se serem 
assinaturas de buracos negros supermassivos. 

As estruturas em rádio do jato de larga es¬ 
cala (Owen et al. 1989; de Gasperin et al. 2012) 
e do núcleo de M87 (Reid et al. 1989; Junor 
et al. 1999; Hada et al. 2016; Mertens et 
al. al. 2016; Kim et ah, 2018b; Walker et ah, 
2018) foram resolvidos em grandes detalhes e em 
mútiplos comprimentos de onda. Além disso, 
a correção do efeito “deslocamento do centro” 
(Blandford & Kónigl 1979), onde a posição apa¬ 


rente do núcleo na frequência de rádio se des¬ 
loca na direção do jato a montante com com¬ 
primento de onda decrescente de maior trans¬ 
parência para auto-absorção síncrotron, indica 
que, em um comprimento de onda de 1,3 mm, 
M87* é coincidente com o buraco negro super- 
massivo (Hada et al. 2011). O envelope do jato 
membro mantém uma forma quase parabólica so¬ 
bre uma vasta gama de distâncias de 10^ Tg até 
~ 20rg (Asada & Nakamura 2012; Hada et al 
2013; Nakamura & Asada 2013; Nakamura et ai. 
2018; Walker et al. 2018). 

Observações de VLBI em 1,3 mm revelam um 
diâmetro da região de emissão de ~ 40/ias, o que 
é comparável à estrutura esperada em escala de 
horizonte (Doeleman et al. 2012; Akiyama et al. 
2015). Essas observações, no entanto, não foram 
capazes de visualizar a sombra do buraco negro 
devido à limitada cobertura da linha de base. 

Com base em três medidas recentes da po¬ 
pulação estelar, adotamos uma distância para 
M87 de 16, 8 ± 0, 8 Mpc (Blakeslee et al. 2009; 
Bird et al. 2010; Cantiello et al. 2018, ver Do¬ 
cumento VI). Usando essa distância e a mode¬ 
lagem do brilho da superfície e da dispersão da 
velocidade das estrelas em comprimentos de onda 
óticos (Gebhardt Sz Thomas 2009; Gebhardt et 
al. 2011), inferimos que a massa de M87* é 
M = 6, 2j'jgg X lO^M© (ver Tabela 9 no Docu¬ 
mento VI). Por outro lado, medições de massa 
modelando a estrutura cinemática do disco de 
gás (Harms et al. 1994; Macchetto et al. 1997) 
produzem M = 3,5j'jg’3 x lO^M© (Walshet al. 
2013, Documento VI). Estas duas estimativas 
de massa, via dinâmica estelar e de gás, pre¬ 
veem um diâmetro teórico da sombra para um 
buraco negro de Schwarzschild de 37,ôl^g/ías e 
21,3l^7/ías, respectivamente. 

3 O Telescópio 


O EHT (Paper H) é um experimento VLBI 
que mede diretamente “visibilidades”, ou com¬ 
ponentes de Eourier, da distribuição de brilho em 
rádio no céu. A medida que a Terra gira, cada par 
de telescópios na rede coleta muitas frequências 
espaciais. O arranjo tem uma resolução angu¬ 
lar nominal de A/L, onde A é o comprimento de 
onda observado e L é o comprimento de linha 
de base máximo projetado entre telescópios no 
arranjo (Thompson et al. 2017). Dessa forma, 
o VLBI cria um telescópio virtual que abrange 
quase todo o diâmetro da Terra. 

Para medir as visibilidades interferométricas, 
os telescópios amplamente separados amostram 



e registram simultânea e coerentemente o campo 
de radiação da fonte. A sincronização usando 
o Sistema de Posicionamento Global (GPS) nor¬ 
malmente atinge o alinhamento em tempo des¬ 
sas gravações em dezenas de nanossegundos. 
Gada estação está equipada com um padrão 
de frequência de maser de hidrogênio. Gom 
as condições atmosféricas durante nossas ob¬ 
servações, o tempo de integração coerente foi ti¬ 
picamente de lOs (ver Figura 2 no Documento 
II). O uso de padrões de frequência de maser de 
hidrogênio em todos os sítios do EHT garante 
coerência em todo o array nesta escala de tempo. 
Após as observações, as gravações são organiza¬ 
das em um local central, alinhadas no tempo, e 
os sinais de cada par de telescópios são correlaci¬ 
onados. 

Embora o VLBI esteja bem estabelecido em 
comprimentos de onda de centímetro e milímetro 
(Boccardi et al. 2017; Thompson et al. 2017) 
e possa ser usado para estudar os ambientes 
próximos de buracos negros (Krichbaum et al. 
1993; Doeleman et al. 2001), a extensão do VLBI 
a um comprimento de onda de 1,3 mm exigiu de¬ 
senvolvimentos técnicos a longo prazo. Desahos 
em comprimentos de onda mais curtos incluem 
aumento do ruído na eletrônica do receptor de 
rádio, mais alta opacidade atmosférica, maiores 
flutuações de fase causadas pela turbulência at¬ 
mosférica e diminuição na ehciência e no tamanho 
dos radiotelescópios nas bandas de observação 
milimétricas e submilimétricas. Iniciado em 2009 
(Doeleman et al. 2009a), o EHT iniciou um pro¬ 
grama para enfrentar esses desahos aumentando 
a sensibilidade do arranjo. O desenvolvimento e 
a implantação de sistemas VLBI de banda larga 
(Whitney et al. 2013; Vertatschitsch et al. 2015) 
levaram a taxas de registro de dados que agora 
excedem as taxas dos arranjos típicos de VLBI 
centimétricas em mais de uma ordem de gran¬ 
deza. Esforços paralelos para apoiar atualizações 
de infraestrutura em locais adicionais de VLBI, 
incluindo o Atacama Large Millimeter / submil- 
limeter Array (ALMA; Matthews et al. 2018; 
Goddi et al. 2019) e o telescópio Atacama Path- 
hnder Experiment (APEX) no Chile (Wagner et 
al. 2015), o Grande Telescópio Milimétrico Al- 
fonso Serrano (LMT) no México (Ortiz-León et 
al. 2016), o telescópio IRAM de 30m no Pico 
Veleta (PV) na Espanha (Greve et al. 1995), o 
Submillimeter Telescope Observatory no Arizona 
(SMT; Baars et al. 1999), o telescópio James 
Clerk Maxwell (JCMT), o Submillimeter Array 
(SMA) no Havaí (Doeleman et al. 2008; Primiani 
et al. 2016; Young et al. 2016), e o Telescópio do 
Polo Sul (SPT) na Antártida (Kim et al. 2018a) 
ampliaram o leque de linhas de base e cobertura 


do EHT e a área geral de coleta do arranjo. Es¬ 
tes desenvolvimentos aumentaram a sensibilidade 
do EHT por um fator de ~ 30 sobre os primeiros 
experimentos que conhrmaram as estruturas em 
escala de horizonte em M87* e Sgr* (Doeleman 
et al. 2008,2012; Akiyama et al. 2015; Johnson 
et al. 2015; Fish et al. 2016; Lu et al. 2018). 

Para as observações em um comprimento de 
onda de 1,3 mm apresentadas aqui, a colaboração 
do EHT utilizou um conjunto global de VLBI de 
oito estações em seis sítios geográhcos. Os com¬ 
primentos de linha de base variaram de 160 a 
10.700 km em direção a M87*, resultando em um 
arranjo com uma resolução teórica de limite de 
difração de ~ 25//as (ver Eiguras 1 e 2 e Docu¬ 
mento H). 

4 Observações, Correlação e Calibração 



Eigura 1: Oito estações da campanha EHT 2017 
em seis localizações geográhcas, conforme visua¬ 
lizadas a partir do plano equatorial. Linhas de 
base cheias representam visibilidade mútua em 
M87* (declinação de -|-12°). As linhas de base 
tracejadas foram usadas para a fonte de cali- 
bração 3C279 (ver Documentos Hl e IV ). 

Observamos M87* com o EHT em 5, 6, 10 e 11 
de abril de 2017. O tempo foi uniformemente de 
bom a excelente, com a mediada da opacidade 
atmosférica noturna no zénite em 230 GHz va¬ 
riando de 0,03 a 0,28 nos diferentes locais. As 
observações foram agendadas com uma série de 
varreduras de três a sete minutos de duração, com 



varreduras em M87* intercaladas com as do qua¬ 
sar 3C 279. O número de varreduras por noite 
obtidas para M87* variou de 7 (10 de abril) a 25 
(6 de abril) como resultado de diferentes horários 
de observação. Uma descrição das observações do 
M87*, suas correlações, calibração e dados finais 
validados é apresentada no Documento III e bre¬ 
vemente resumida aqui. 

Em cada estação, o sinal astronômico em 
ambas as polarizações e em duas bandas de 
frequências largas adjacentes de 2 GHz centra¬ 
das em 227,1 e 229,1 GHz foram convertidas para 
banda base usando técnicas heteródinas padrão, 
depois digitalizadas e gravadas a uma taxa to¬ 
tal de 32 Gbps. A correlação dos dados foi re¬ 
alizada usando um correlacionador de software 
(Deller et al. 2007) no MIT Haystack Observa- 
tory e no Max-Planck-Institut für Radioastrono- 
mie, cada um lidando com uma das duas bandas 
de frequência. As diferenças entre os dois cor- 
relacionadores independentes mostraram-se in¬ 
significantes através da troca de algumas varre¬ 
duras idênticas para comparação cruzada. Na 
correlação, os sinais foram alinhados a uma re¬ 
ferência de tempo comum usando uma geometria 
da Terra a priori e um modelo de relógio. 

Uma etapa subsequente de ajuste de frente 
de onda identificou detecções na potência de 
sinal correlacionada quando calibrando a fase 
dos dados para atrasos residuais e efeitos at¬ 
mosféricos. A utilização do ALMA como 
estação de referência altamente sensível permitiu 
correções críticas para distorções ionosféricas e 
troposféricas nos outros locais. O ajuste de frente 
de onda foi realizado com três pipelines automa¬ 
tizadas independentes, cada uma adaptada às ca¬ 
racterísticas específicas das observações do EHT, 
tais como largura de banda larga, suscetibilidade 
à turbulência atmosférica e heterogeneidade do 
arranjo (Blackburn et al., 2019; Janssen et al., 
2019, Paper III). Os pipelines fizeram uso de soft¬ 
ware padrão para o processamento de dados de 
rádio-interferometria (Greisen 2003; Whitney et 
al. 2004; McMullin et al. 2007, IM van Bemmel 
et al. 2019, em preparação). 

Os dados das pipelines de ajuste de frente de 
onda foram escalonados dos coeficientes de cor¬ 
relação para uma escala de densidade de fluxo 
físico uniforme (em Jansky) usando uma esti¬ 
mativa inicial independente da sensibilidade de 
cada telescópio. As precisões das sensibilidades 
da estação foram estimadas entre 5% a 10% em 
amplitude, embora certas perdas não caracteriza¬ 
das (por exemplo, de mau apontamento ou foco) 
possam exceder a estimativa de erro. Assumindo 
valores de densidade de fluxo total derivados de 



Eigura 2: Acima: Cobertura {u, v) para M87*, 
agregada ao longo de todos os quatro dias das 
observações. As coordenadas (rt, v) para cada 
par de antenas são o comprimento da linha de 
base projectada na fonte em unidades do com¬ 
primento de onda de observação A e são dadas 
para pares conjugados. Linhas de base para o 
ALMA/APEX e para o JCMT/SMA são redun¬ 
dantes. Linhas circulares pontilhadas indicam 
comprimentos de linha de base correspondentes 
a espaçamentos de franja de 50 e 25^as. In¬ 
ferior: amplitudes de visibilidade calibradas fi¬ 
nais de M87* em função do comprimento de li¬ 
nha de base projetado em II de abril. As linhas 
de base redundantes para APEX e JCMT são 
plotadas como losangos. Barras de erro corres¬ 
pondem a incertezas térmicas (estatísticas). A 
transformada de Fourier de um modelo de anel 
fino simétrico azimutalmente com diâmetro de 
46^as é também mostrado com uma linha tra¬ 
cejada para comparação. 














dados interferométricos do ALMA (Goddi et al. 
2019) e utilizando redundância de matriz via cali- 
bração de rede (Paper III), refinamos a calibração 
de amplitude absoluta de telescópios que estão 
posicionados e possuem linhas de base redundan¬ 
tes, ou seja, ALMA/APEX e JCMT/SMA. 

A mediana da razão sinal/ruído da média da 
varredura para M87* foi > 10 em linhas de base 
não-ALMA e > 100 em linhas de base para o 
ALMA, levando a pequenos erros estatísticos na 
amplitude e fase de visibilidade. As comparações 
entre as três pipelines independentes, as duas po¬ 
larizações e as duas bandas de frequência permi¬ 
tiram a estimativa de erros de linha de base sis¬ 
temáticos de cerca de 1° na fase de visibilidade e 
de 2% nas amplitudes de visibilidade. Esses pe¬ 
quenos erros limitantes permanecem mesmo de¬ 
pois do ajuste das sensibilidades das estações e 
fases desconhecidas da estação através da autoca- 
libração (Pearson Sz Readhead 1984) e afetam as 
quantidades de fechamento interferométrico (Ro- 
gers et al., 1974; Readhead et al., 1980). Após 
a validação de dados e comparações de pipeline, 
uma única saída da pipeline foi designada como o 
conjunto de dados primário da primeira liberação 
de dados científicos do EHT e usada para resul¬ 
tados subsequentes, enquanto as saídas dos ou¬ 
tros dois pipelines ofereceram suporte a conjun¬ 
tos de dados de validação. As visibilidades com¬ 
plexas hnais calibradas V{u,v) correspondem as 
componentes de Pourier da distribuição de bri¬ 
lho no céu na frequência espacial {u, v) determi¬ 
nada pela linha de base projetada expressa em 
unidades do comprimento de onda de observação 
(van Cittert 1934; Thompson et al. 2017). A Fi¬ 
gura 2 mostra a cobertura (u, v) e as amplitudes 
de visibilidades calibradas de M87* para 11 de 
abril. As amplitudes de visibilidade se asseme¬ 
lham às de um anel hno (isto é, uma função de 
Bessel Jq; veja a Figura 10.12 em Thompson et 
al. 2017). Esse modelo de anel com diâmetro de 
46 fias tem um primeiro nulo a 3,4 GA, combi¬ 
nando o mínimo na densidade de fluxo observada 
e é consistente com uma densidade de fluxo re¬ 
duzida na linha de base mais longa do Hawaih- 
Espanha (JGMT / SMA-PV) prõximo de 8 GA. 
Este modelo particular de anel, mostrado com 
uma linha tracejada no painel inferior da Figura 
2, é apenas ilustrativo e não se ajusta a todos 
os recursos nos dados. Primeiro, as amplitudes 
de visibilidade nas linhas de base de VLBI mais 
curtas sugerem que cerca de metade da densidade 
de fluxo compacta vista na linha de base ALMA- 
APEX de ~2 km é resolvida pelo feixe interfe¬ 
rométrico (Paper IV). Em segundo lugar, dife¬ 
renças na profundidade do primeiro mínimo em 
função da orientação, bem como fases de fecha¬ 


mento medidas altamente diferentes de zero, in¬ 
dicam algum grau de assimetria na fonte (Papers 
III , VI). Finalmente, as amplitudes de visibili¬ 
dade representam apenas metade da informação 
disponível para nós. Em seguida, exploraremos 
imagens e modelos geométricos mais complexos 
que podem se ajustar às amplitudes e fases de 
visibilidade medidas. 


5 Imagens e Recursos 


Nós reconstruímos imagens calibrando as visi¬ 
bilidades do EHT, que fornecem resultados inde¬ 
pendentes dos modelos (Papel IV ). No entanto, 
existem dois grandes desahos na reconstrução de 
imagens a partir de dados EHT. Em primeiro 
lugar, as linhas de base EHT amostram uma 
gama limitada de frequências espaciais, corres¬ 
pondendo a escalas angulares entre 25 e 160 fias. 
Gomo o plano (n, v) é escassamente amostrado 
(Figura 2), o problema inverso é sub- vinculado. 
Segundo, as visibilidades medidas não possuem 
calibração de fase absoluta e podem ter incerte¬ 
zas de calibração de grande amplitude. 

Para lidar com esses desahos, os algoritmos 
de imagens incorporam suposiçóes e vínculos adi¬ 
cionais que são projetadas para produzir ima¬ 
gens que são hsicamente plausíveis (por exem¬ 
plo, positivas e compactas) ou conservadoras (por 
exemplo, suaves), embora permaneçam consis¬ 
tentes com os dados. Nós exploramos duas clas¬ 
ses de algoritmos para reconstruir imagens a par¬ 
tir de dados do EHT. A primeira classe de algo¬ 
ritmos é a abordagem tradicional GLEAN usada 
na rádio interferometrica (por exemplo, Hõgbom 
1974; Glark 1980). GLEAN é uma abordagem 
de modelagem inversa que faz a convolução in¬ 
versa da função de espalhamento pontual a par¬ 
tir das visibilidades no espaço das transforma¬ 
das de Fourier. Ao aplicar GLEAN, é necessário 
auto-calibrar iterativamente os dados entre as 
rodadas de geração de imagens para resolver 
erros de fase e amplitude variáveis no tempo 
nos dados. A segunda classe de algoritmos é a 
chamada Máxima Verossimilhança Regularizada 
(RML na sigla em inglês; por exemplo, Narayan e 
Nityananda, 1986; Wiaux et al., 2009; Thiébaut, 
2013). A RML é uma abordagem de modela¬ 
gem avançada que procura uma imagem que não 
seja apenas consistente com os dados observa¬ 
dos, mas também favoreça propriedades de ima¬ 
gem especíhcas (por exemplo, suavidade ou com¬ 
pactação). Assim como com GLEAN, os métodos 
RML normalmente iteram entre a geração de 
imagens e a autocalibração, embora também pos¬ 
sam ser usados para imagens diretamente em 



quantidades robustas de fechamento imunes a 
erros de calibração baseados em estações. Os 
métodos RML foram extensivamente desenvol¬ 
vidos para o EHT (por exemplo, Honma et al. 
2014; Bouman et al. 2016; Akiyama et al. 2017; 
Chael et al. 2018b; ver também Paper IV). 

Todo algoritmo de imagem possui uma varie¬ 
dade de parâmetros livres que podem afetar sig¬ 
nificativamente a imagem final. Adotamos uma 
abordagem de imagem de dois estágios para con¬ 
trolar e avaliar vieses nas reconstruções a par¬ 
tir de nossas escolhas desses parâmetros. No 
primeiro estágio, quatro equipes trabalharam de 
forma independente para reconstruir as primei¬ 
ras imagens do EHT do M87* usando uma li¬ 
beração antecipada de dados de engenharia. As 
equipes trabalharam sem interação para minimi¬ 
zar o viés compartilhado, mas cada uma delas 
produziu uma imagem com uma característica 
proeminente semelhante: um anel de diâmetro 
~ 38 —44//as, com brilho relevante na direção sul 
(veja a Eigura 4 no Paper IV). 

No segundo estágio de imagem, desenvolve¬ 
mos três pipelines de imagem, cada um usando 
um pacote de software diferente e uma metodo¬ 
logia associada. Cada pipeline pesquisou uma 
série de parâmetros de imagem, produzindo ~ 
10^ a 10^ imagens de diferentes combinações de 
parâmetros. Determinamos um conjunto de me¬ 
lhores combinações de parâmetros que tanto pro¬ 
duziram imagens de M87* que eram consisten¬ 
tes com os dados observados e que reconstruíram 
imagens precisas de conjuntos de dados sintéticos 
correspondentes a quatro modelos geométricos 
conhecidos (anel, crescente, disco preenchido e 
fonte dupla assimétrica). Para todos os pipelines, 
os melhores conjuntos de dados mostraram um 
anel assimétrico com um diâmetro de ~ 40/ias, 
com diferenças surgindo principalmente nas re¬ 
soluções angulares efetivas obtidas por diferentes 
métodos. 

Para cada pipeline, determinamos a com¬ 
binação única de parâmetros de imagem fidu- 
cial fora dos melhores conjuntos que tiveram me¬ 
lhor desempenho em todos os conjuntos de dados 
sintéticos e para cada metodologia de imagem as¬ 
sociada (ver Figura 11 no Paper IV). Como as re¬ 
soluções angulares das imagens reconstruídas va¬ 
riam entre os pipelines, nós desfocamos cada ima¬ 
gem com uma Gaussiana circular para uma re¬ 
solução angular conservadora comum de 20 fias. 
A parte superior da Figura 3 mostra uma ima¬ 
gem de M87* em 11 de abril obtida pela média 
das imagens fiduciais desfocadas das três pipeli¬ 
nes. A imagem é dominada por um anel com um 
perfil azimutal assimétrico que é orientado em 


um ângulo de posição ~ 170° na direção leste, 
partindo do Norte. Embora o ângulo de posição 
medido aumente ~ 20° entre os dois primeiros 
dias e os dois últimos dias, as características da 
imagem são amplamente consistentes nos diferen¬ 
tes métodos de imagem e nos quatro dias de ob¬ 
servação. Isso é mostrado na parte inferior da 
Figura 3, que relata as imagens em dias diferen¬ 
tes (veja também a Figura 15 no Documento IV). 
Estes resultados são também consistentes com os 
obtidos a partir do ajuste do domínio de visi¬ 
bilidade dos modelos de magneto-hidrodinâmica 
geométrica e relatividade geral ou GRMHD (Pa¬ 
per VI). 

6 Modelagem Teórica 


A aparência do M87* foi modelada com su¬ 
cesso usando simulações GRMHD, que descre¬ 
vem um disco magnetizado quente e turbulento 
orbitando um buraco negro de Kerr. Eles na¬ 
turalmente produzem um jato potente e podem 
explicar a distribuição de energia espectral de 
banda larga observada nos LLAGNs. Em um 
comprimento de onda de 1,3 mm, e como obser¬ 
vado aqui, as simulações também prevêem uma 
sombra e um anel de emissão assimétrica. Este 
último não coincide necessariamente com a órbita 
circular estável mais interna, ou ISGO, e está 
relacionado ao anel de fótons com lente. Para 
explorar este cenário em grande detalhe, nós 
construímos uma biblioteca de imagens sintéticas 
(Biblioteca de Imagens) descrevendo fluxos de 
acreção magnetizada em buracos negros em RG* 
(Paper V). As imagens em si são produzidas 
a partir de uma biblioteca de simulações (Bi¬ 
blioteca de Simulações) coletando os resultados 
de quatro códigos resolvendo as equações do 
GRMHD (Gammie et al. 2003; Sadowski et al. 
2014; Porth et al. 2017; Liska et al. 2018). Os 
elementos da Biblioteca de Simulações foram aco¬ 
plados a três traços de raios-relativísticos gerais 
e diferentes códigos de transferência radioativa 
(GRRT, Bronzwaer et al. 2018; Moscibrodzka 
&: Gammie 2018; Z. Younsi et al. 2019, em pre¬ 
paração). Nós nos limitamos a fornecer aqui uma 
breve descrição das configurações iniciais e dos 
cenários físicos explorados nas simulações; veja o 
artigo V para detalhes sobre os códigos GRMHD 
e GRRT, que foram co-validados em precisão e 
consistência (Gold et al. 2019; Porth et al. 2019). 

Uma típica simulação GRMHD na biblioteca 
é caracterizada por dois parâmetros: o spin adi- 

* Espaços espaciais mais exóticos, como buracos negros 
de dilaton, estrelas de bósons e gravastares, também foram 
considerados (Artigo V). 




mensional a* = Jc/GM"^, onde J e M são, res¬ 
pectivamente, o momento angular de rotação e 
a massa do buraco negro, e o fluxo magnético 
adimensional líquido sobre o horizonte de evento 

(j) = ^/{M, onde $ e M são o fluxo 
magnético e o fluxo de massa (ou taxa de acreção) 
ao longo do horizonte, respectivamente. Como as 
simulações do GRMHD escalam com a massa do 
buraco negro, M é definido apenas no momento 
de produzir as imagens sintéticas com os códigos 
GRRT. O fluxo magnético é geralmente diferente 
de zero porque o campo magnético fica preso no 
buraco negro pelo fluxo de acreção e é sustentado 
pelas correntes no plasma circundante. 

Esses dois parâmetros nos permitem descre¬ 
ver discos de acreção que são prógrados (a* > 
0) ou retrógrados (a* < 0) em relação ao 

eixo de rotação do buraco negro, e cujos flu¬ 
xos de acréscimo são “SANE” (de “Standard and 
Evolução Normal”, Narayan et al. 2012) com 
(?!) ~ 1, ou “MAD” (de “Disco Preso Magnetica¬ 
mente”, Narayan et al. 2003) com cj) 15^. Em 
essência, os fluxos de acreção SANE são carac¬ 
terizados por um fluxo magnético adimensional 
moderado e resultam de campos magnéticos ini¬ 
ciais menores que os dos fluxos MAD. Além disso, 
os ângulos de abertura do funil magnético nos 
fluxos SANE são geralmente menores que os dos 
fluxos MAD. Variando a* e (j), realizamos 43 lon¬ 
gas simulações tridimensionais de alta resolução, 
cobrindo bem as propriedades físicas dos fluxos 
de acreção magnetizada nos buracos negros de 
Kerr. 

Todas as simulações GRMHD foram inicia- 
lizadas com um toro fracamente magnetizados 
em órbita de um buraco negro e conduzido para 
um estado de turbulência por instabilidades, in¬ 
cluindo a instabilidade magnetorotacional (Bal- 
bus & Hawley 1991), rapidamente atingindo um 
estado quase estacionãrio. Uma vez concluída 
a simulação, as propriedades de fluxo relevantes 
em diferentes momentos são coletadas para serem 
empregadas para o pós-processamento adicional 
dos códigos da GRRT. A geração de imagens de 
síntese requer, além das propriedades do fluido 
(de campo magnético, do campo de velocidade e 
densidade de massa de repouso), também os co¬ 
eficientes de emissão e absorção, a inclinação i 
(o ângulo entre o vetor de momento angular de 
fluxo de acreção e a linha de visão), o ângulo 
de posição PA (o ângulo a partir do Norte na 
direção leste, isto é, no sentido anti-horário nas 
nossas imagens, da projeção no céu do momento 
angular do fluxo de acreção), a massa do buraco 

^Nós aqui usamos unidades Heaviside, onde um fator 
de -\/47r é absorvido no definição do campo magnético. 
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Figura 3: Acima: Imagem EHT do M87* das ob¬ 
servações de 11 de abril de 2017 como um exem¬ 
plo representativo das imagens coletadas na cam¬ 
panha de 2017. A imagem é a média de três di¬ 
ferentes métodos de imageamento após convoluir 
cada uma com um kernel Gaussiano circular para 
fornecer resoluções correspondentes. O maior dos 
três kernels (20//as de largura a meia-altura) é 
mostrado no canto inferior direito. A imagem é 
mostrada em unidades de temperatura de brilho 
Tb = 5A^/2 A:bD, onde S é a densidade de fluxo, 
A é o comprimento de onda observado, ks é a 
constante de Boltzmann e D é o ângulo sólido do 
elemento de resolução. Abaixo: imagens seme¬ 
lhantes tiradas em diferentes dias, mostrando a 
estabilidade da estrutura básica da imagem e as 
semelhanças nos diferentes dias. Norte acima e 
leste à esquerda. 

negro Mea distância D ao observador. 

Porque se acredita que os fótons no compri¬ 
mento de onda de 1,3 mm observados pelo EHT 
sejam produzidos por emissão síncrotron, cujos 
coeficientes de absorção e emissão dependem da 
função de distribuição de elétrons, consideramos 
o plasma como composto de elétrons e íons que 










têm a mesma temperatura nas regiões magne¬ 
ticamente dominadas do fluxo (funil), mas têm 
uma temperatura substancialmente diferente nas 
regiões dominadas por gás (plano intermediário 
do disco). Em particular, consideramos o plasma 
como composto de íons não relativísticos com 
temperatura Tj e elétrons relativísticos com tem¬ 
peratura Te- Uma prescrição simples para a razão 
das temperaturas das duas espécies pode então 
ser imposta em termos de um único parâmetro 
(Moscibrodzka et al. 2016), tal que a maior parte 
da emissão provém de regiões magnetizadas fra¬ 
camente (pequenas Raito, Te ~ Ti/Raito) ou for- 
temente magnetizadas (grandes Raito , Tg Ti). 
Nos modelos SANE, o disco (jato) é magneti¬ 
camente fraco (forte), de modo que os modelos 
de baixa (alta) Raito produzem a maior parte da 
emissão no disco (jato). Nos modelos MAD, exis¬ 
tem regiões fortemente magnetizadas em todos os 
lugares e a emissão é principalmente do meio do 
disco. Embora esta prescrição não seja a única 
possível, ela tem a vantagem de ser simples, su¬ 
ficientemente genérica e robusta (ver Documento 
V para uma discussão de partículas não térmicas 
e resfriamento radiativo). 


Tabela 1 

Parâmetros de M87* 


Parâmetro 

Estimativa 

Diâmetro do anel“ d 

42 ± 3/ras 

Largura do anel “ 

< 20/ias 

Contraste crescente*’ 

10:1 

Razão axial “ 

4:1 

Angulo de posição 

150° - 200° E do N 

Og = GM/Dc^ 

3,8 ± 0,4/ias 

a = d/6g 

ii+O.S 

M 

6, 5 ±0,7 X 10®M© 

Parâmetro 

Estimativa prévia 


(16, 8 ±0,8) Mpc 

M(estrelas)® 

6.2lJ;i X 10®M© 

M(gases)® 

S.Slg;® X 10®M© 


Notes. 

“ Derivado do domínio da imagem. 

^ Derivado do ajuste do modelo do crescente. 

A massa e os erros sistemáticos são médias dos três 
métodos (modelos geométricos, modelos GRMHD e 
extração do anel de domínio de imagem). 

^ O valor exato depende do método usado para 
extrair d, que é refletido no intervalo fornecido. 

® Rederivado das distribuições de probabilidade 
(Paper VI). 


Como cada simulação GRMHD pode ser 
usada para descrever vários cenários físicos di¬ 
ferentes, alterando a função de distribuição de 


elétrons prescrita, usamos a Biblioteca de Si¬ 
mulação para gerar mais de 420 cenários físicos 
diferentes. Cada cenário é então usado para gerar 
centenas de instantâneos em diferentes momentos 
da simulação, levando a mais de 62.000 objetos 
na Biblioteca de Imagens. A partir das imagens, 
criamos observáveis do modelo que correspondem 
ao cronograma de observação do EHT e as com¬ 
paramos com as observáveis de VLBI medidas, 
conforme detalhado no Documento VI. 

Como exemplo, a linha superior da Figura 4 
mostra três instantâneos do modelo GRMHD da 
Biblioteca de Imagens com diferentes rotações e 
tipo de fluxo, que melhor ajustam as fases de 
fechamento e as amplitudes dos dados de 11 de 
abril. Para esses modelos, produzimos visibilida¬ 
des simuladas para os parâmetros do cronograma 
e do clima de 11 de abril e as calibramos com 
um pipeline sintético de geração de dados e cali- 
bração (Blecher et al. 2017; Janssen et al. 2019; 
Roelofs et ah, 2019a). Os dados simulados produ¬ 
ziram então imagens com o mesmo pipeline usado 
para as imagens observadas. As semelhanças en¬ 
tre as imagens simuladas (linha inferior da Fi¬ 
gura 4 ) e as imagens observadas (Figura 3) são 
notáveis. 

No geral, ao combinar todas as informações 
contidas na Biblioteca de Simulação e na Bi¬ 
blioteca de Imagens, a origem física das carac¬ 
terísticas de emissão da imagem observada em 
M87* pode ser resumida da forma que segue. 

Primeiro, a imagem observada é consistente 
com a hipótese de que é produzida por um fluxo 
de acreção magnetizado orbitando dentro de al¬ 
guns rg do horizonte de eventos de um buraco 
negro de Kerr. O anel assimétrico é produzido 
por uma combinação de lentes gravitacionais for¬ 
tes e feixes relativísticos, enquanto a depressão 
do fluxo central é a assinatura observacional da 
sombra do buraco negro. Curiosamente, todos os 
modelos de acreção são consistentes com a ima¬ 
gem EHT, exceto para os modelos MAD com 
a* = —0.94, que não conseguem produzir ima¬ 
gens suhcientemente estáveis (isto é, a variação 
entre instantâneos é demasiado grande para ser 
estatisticamente consistente com o imagem ob¬ 
servada) . 

Em segundo lugar, a assimetria norte-sul no 
anel de emissão é controlada pela rotação do bu¬ 
raco negro e pode ser usada para deduzir sua ori¬ 
entação. Nos modelos de disco co-rotante (onde 
o momento angular da matéria e do buraco negro 
estão alinhados), a parede do funil, ou bainha do 
jato, gira com o disco e o buraco negro; nos mo¬ 
delos de disco com rotação contrária, a parede do 
funil luminoso gira com o buraco negro, mas con- 



tra o disco. A assimetria norte-sul é consistente 
com modelos em que o giro do buraco negro está 
apontando para longe da Terra e inconsistente 
com os modelos em que o vetor de giro aponta 
para a Terra. 

Terceiro, adotando uma inclinação de 17° en¬ 
tre o jato que se aproxima e a linha de visada 
(Walker et al. 2018 ), a orientação oeste do jato e 
um modelo de disco corotante, a matéria na parte 
inferior da imagem está se movendo em direção 
ao observador (rotação no sentido horário visto 
da Terra). Isto é consistente com a rotação do 
gás ionizado em escalas de 20 pc, ou seja, TOOOr^ 
(Ford et al. 1994; Walsh et al. 2013 ) e com o 
senso inferido de rotação de observações VLBI a 
7 mm (Walker et al. al. 2018 ). 

Finalmente, modelos com um a* = 0 são des¬ 
favorecidos pela exigência observacional muito 
conservadora de que a potência do jato seja 
Pjet > 10^^ erg s“^. Além disso, nos modelos que 
produzem um jato suficientemente potente, ele é 
alimentado pela extração da energia de rotação 
do buraco negro através de mecanismos seme¬ 
lhantes ao processo de Blandford-Znajek. 
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Figura 4: Parte superior: três exemplos de 
modelos de alguns dos instantâneos de melhor 
ajuste da biblioteca de imagens de simulações 
do GRMHD para o dia 11 de abril, correspon¬ 
dentes a diferentes parâmetros de giro e fluxos 
de acreção. Parte inferior: os mesmos mode¬ 
los teóricos, processados por meio de um pipe- 
line de simulação VLBI com as mesmas datas, 
características do telescópio e parâmetros meteo¬ 
rológicos, como na edição de 11 de abril, e image- 
ados da mesma forma que a Figura 3. Note que, 
embora o ajuste às observações seja igualmente 
bom nos três casos, eles se referem a cenários 
físicos radicalmente diferentes; isso destaca que 
um único bom ajuste não implica que um mo¬ 
delo seja preferido em detrimento de outros (veja 
o Documento V). 


7 Comparação de Modelos e Estimativa 
de Parâmetros 


No Documento VI, a massa do buraco negro é 
derivada dos ajustes aos dados visíveis dos mode¬ 
los geométricos e GRMHD, bem como das medi¬ 
das do diâmetro do anel no domínio da imagem. 
Nossas medições permanecem consistentes entre 
metodologias, algoritmos, representações de da¬ 
dos e conjuntos de dados observados. 

Motivado pelas estruturas anelares de emissão 
assimétricas vistas nas imagens reconstruídas 
(Seção 5) e estruturas de emissão similares vis¬ 
tas nas imagens das simulações de GRMHD 
(Seção 6), desenvolvemos uma família de mode¬ 
los geométricos em crescente (ver, por exemplo, 
Kamruddin & Dexter 2013) para comparar dire¬ 
tamente com os dados visíveis. Utilizamos dois 
algoritmos de inferência Bayesiana e demonstra¬ 
mos que tais modelos crescentes são estatistica¬ 
mente preferidos em comparação com outros mo¬ 
delos geométricos comparativamente complexos 
que exploramos. Descobrimos que os modelos 
crescentes fornecem ajustes aos dados que são 
estatisticamente comparáveis aos das imagens re¬ 
construídas apresentadas na Seção 5, permitindo- 
nos determinar os parâmetros básicos dos cres¬ 
centes. Os modelos de melhor ajuste para o 
anel de emissão assimétrica têm diâmetros de 
43 ± 0,9/uas e larguras fracionais em relação ao 
diâmetro de <0,5. A emissão cai abruptamente 
no interior do anel de emissão assimétrica com a 
depressão central tendo um brilho <10% do bri¬ 
lho médio no anel. 

Os diâmetros dos modelos crescentes 
geométricos medem os tamanhos característicos 
das regiões emissoras que circundam as sombras 
e não os tamanhos das próprias sombras (ver, 
por exemplo, Psaltis et al. 2015; Johannsen et 
al. 2016; Kuramochi et al. 2018, para potenciais 
vieses). 

Nós modelamos o diâmetro angular cescente 
d em termos do raio e distância gravitacional 
9g = GM/c^D, como d = a9g, onde a é uma 
função da rotação, inclinação, e Raito (o — 9.6 — 
10.4) corresponde a emissão a partir apenas do 
anel de fotons lenteado. Nós calibramos a ajus¬ 
tando os modelos crescentes geométricos a um 
grande número de dados observáveis gerados a 
partir da Biblioteca de Imagens. Também po¬ 
demos ajustar os observáveis modelados gerados 
a partir da Biblioteca de Imagens para os da¬ 
dos M87*, o que nos permite medir 9g direta¬ 
mente. No entanto, tal procedimento é com¬ 
plicado pela natureza estocástica da emissão no 
fluxo de acreção (ver, por exemplo, Kim et al. 











2016). Para dar conta dessa estrutura turbu¬ 
lenta, desenvolvemos um formalismo e vários al¬ 
goritmos que estimam as estatísticas dos com¬ 
ponentes estocásticos usando conjuntos de ima¬ 
gens de simulações individuais de GRMHD. Des¬ 
cobrimos que os dados observáveis não são in¬ 
consistentes com as simulações do GRMHD. Con¬ 
cluímos que os parâmetros do modelo recuperado 
são consistentes entre os algoritmos. 

Finalmente, extraímos o diâmetro, a largura 
e a forma do anel diretamente das imagens re¬ 
construídas (consulte a Seção 5). Os resul¬ 
tados são consistentes com as estimativas de 
parâmetros de modelos crescentes geométricos. 
Seguindo o mesmo procedimento de calibração 
GRMHD, infere-se valores de 6g e a de probabi¬ 
lidade mãxima regulares e imagens reconstruídas 
CLEAN. 

Combinando resultados de todos os métodos, 
medimos diâmetros de região de emissão de 42 ± 
3/uas, tamanhos angulares do raio gravitacional 
6g = 3,8 ± 0,4^as, e fatores de escala no in¬ 
tervalo alpha = 10,7-11,5, com erros associa¬ 
dos de ~ 10%. Para a distância de 16, 8 ± 0,8 
Mpc adotada aqui, a massa do buraco negro é 
M = (6,5 ± 0,7) X 1O®M0; o erro sistemático 
refere-se ao nível de confiança de 68% e é muito 
maior do que o erro estatístico de 0, 2 x 1O®M0. 
Além disso, ao traçar o pico da emissão no anel, 
podemos determinar a forma da imagem e ob¬ 
ter uma relação entre o eixo principal e o menor 
do anel, que é aproximadamente menor que 4:3; 
isto corresponde a um desvio aproximadamente 
menor que 10% da circularidade em termos de 
distância rms a partir de um raio médio. 

A Tabela 1 resume os parâmetros medidos ca¬ 
racterísticos da imagem e as propriedades de bu¬ 
raco negro inferidas com base nos dados de todas 
as bandas e todos os dias combinados. A massa 
do buraco negro inferida favorece fortemente a 
medição baseada na dinâmica estelar (Gebhardt 
et al. 2011). O tamanho, a assimetria, o con¬ 
traste de brilho e a circularidade das imagens 
reconstruídas e dos modelos geométricos, bem 
como o sucesso das simulações GRMHD na des¬ 
crição dos dados interferométricos, são consisten¬ 
tes com as imagens EHT do M87* sendo associa¬ 
das à emissão fortemente lenteada da vizinhança 
de um buraco negro de Kerr. 

8 Discussão 


Vários elementos reforçam a robustez de nossa 
imagem e a conclusão de que ela é consistente 
com a sombra de um buraco negro como pre¬ 
dito pela RG. Primeiro, nossa análise usou várias 


técnicas independentes de calibração e imagem, 
bem como quatro conjuntos de dados indepen¬ 
dentes realizados em quatro dias diferentes em 
duas bandas de frequência separadas. Em se¬ 
gundo lugar, a estrutura da imagem corresponde 
bem às previsões anteriores (Dexter et al. 2012; 
Moscibrodzka et al. 2016) e é bem reproduzida 
pelo nosso esforço de modelagem extensiva apre¬ 
sentado na Seção 6. Terceiro, porque a nossa 
medição da massa do buraco negro em M87* não 
é inconsistente com todas as medições anterio¬ 
res de massa, e isso nos permite concluir que a 
hipótese nula da métrica de Kerr (Psaltis et al. 
2015; Johannsen et al. 2016), ou seja, a suposição 
de que o buraco negro é descrito pela métrica 
de Kerr, não foi violada. Quarto, o anel de 
emissão observado reconstruído em nossas ima¬ 
gens é próximo ao circular com uma razão axial 
aproximadamente < 4:3; da mesma forma, as 
imagens médias de nossas simulações GRMHD 
também mostram uma forma circular. Depois de 
associar à forma da sombra um desvio da circu¬ 
laridade - medido em termos de distância média 
quadrática de um raio médio na imagem - que é 
aproximadamente menor que 10%, podemos defi¬ 
nir um limite inicial de ordem quatro em desvios 
relativos do momento de quadrupólo do valor de 
Kerr (Johannsen &: Psaltis 2010). Dito isto de 
forma diferente, se Q é o momento quadrupolo 
de um buraco negro de Kerr e /S.Q o desvio como 
deduzido a partir da circularidade, nossa medida 
e o fato de que o ângulo de inclinação é assumido 
como sendo pequeno implica que /S.Q/Q < 4 
{/S.Q/Q = e em Johannsen & Psaltis 2010). 

Por fim, ao comparar as amplitudes de visi¬ 
bilidade do M87* com os dados de 2009 e 2012 
(Doeleman et al. 2012; Akiyama et al. 2015), o 
tamanho total do núcleo de rádio em um compri¬ 
mento de onda de 1,3 mm não mudou significa¬ 
tivamente, apesar da variabilidade na densidade 
total de fluxo. Essa estabilidade é consistente 
com a expectativa de que o tamanho da sombra 
é uma característica ligada à massa do buraco ne¬ 
gro e não às propriedades de um fluxo de plasma 
variãvel. 

Também é fácil rejeitar algumas inter¬ 
pretações astrofísicas alternativas. Por exemplo, 
é improvável que a imagem seja produzida por 
uma característica de jato, já que observações de 
várias épocas de VLBI do jato de plasma em M87 
(Walker et al. 2018) em escalas fora do hori¬ 
zonte não mostram anéis circulares. O mesmo 
é tipicamente verdadeiro para jatos AGN em 
grandes pesquisas de VLBI (Lister et al. 2018). 
Similarmente, se o anel aparente fosse um ali¬ 
nhamento aleatório de bolhas de emissão, eles 



também deveriam ter se deslocado em velocida¬ 
des relativísticas, isto é, em ~ òfias dia (Kim 
et al. 2018b), levando a mudanças na estru¬ 
tura e nos tamanhos mensuráveis. Os modelos 
GRMHD de cones de jato oco podem mostrar, 
sob condições extremas, características de anel 
estáveis (Pu et al. 2017), mas este efeito está in¬ 
cluído em certa medida em nossa Biblioteca de 
Simulação para modelos com Raito > 10. Final¬ 
mente, um anel de Einstein formado por lente 
gravitacional de uma região brilhante no contra- 
jato exigiria um alinhamento bem ajustado e um 
tamanho maior do que o medido em 2012 e 2009. 

Ao mesmo tempo, é mais difícil descartar al¬ 
ternativas aos buracos negros na RG, porque uma 
sombra pode ser produzida por qualquer objeto 
compacto com um espaço-tempo caracterizado 
por órbitas de fótons circulares instáveis (Mizuno 
et al., 2018). De fato, embora a métrica de Kerr 
continue a ser uma solução em algumas teorias 
alternativas da gravidade (Barausse & Sotiriou 
2008; Psaltis et al. 2008), existem soluções de 
buracos negros não-Kerr em uma variedade des¬ 
sas teorias modificadas (Berti et al. 2015). Além 
disso, alternativas exóticas aos buracos negros, 
como singularidades nuas (Shaikh et al. 2019), 
estrelas do bóson (Kaup 1968; Liebling k, Pa- 
lenzuela 2012) e gravastars (Mazur k Mottola 
2004; Chirenti k Rezzolla 2007), são soluções 
admissíveis dentro da RG e fornecem modelos 
concretos, ainda que controversos. Alguns desses 
objetos compactos exóticos já podem ser mostra¬ 
dos como incompatíveis com nossas observações, 
dada a nossa massa máxima anterior. Por exem¬ 
plo, as sombras de singularidades nuas associadas 
aos espaços-tempos de Kerr |a*| > 1 são substan¬ 
cialmente menores e muito assimétricas compara¬ 
das com as dos buracos negros de Kerr (Bambi 
k Freese 2009). Além disso, alguns tipos comuns 
de buracos de minhoca (Bambi 2013) prevêem 
sombras muito menores do que as que medimos. 

No entanto, outros candidatos a objeto com¬ 
pacto precisam ser analisados com mais cuidado. 
As estrelas de Boson são um exemplo de obje¬ 
tos compactos com órbitas circulares de fótons, 
mas sem uma superfície ou um horizonte de even¬ 
tos. Em tal espaço-tempo, geodésicas nulas são 
redirecionadas para fora em direção a observado¬ 
res distantes (Cunha et al. 2016), de modo que 
a sombra pode em princípio ser preenchida com 
emissão de imagens de fontes de rádio distantes 
gerando uma imagem espelhada complexa do céu. 
Mais importante ainda, os fluxos de acreção nas 
estrelas do bóson se comportam de maneira dife¬ 
rente, pois não produzem jatos, mas protelam os 
toros de acreção que os distinguem dos buracos 


negros (Olivares et ah, 2019). Gravastars forne¬ 
cem exemplos de objetos compactos com órbitas 
instáveis de fótons e uma superfície rígida, mas 
não um horizonte de eventos. Neste caso, en¬ 
quanto uma única imagem do fluxo de acreção 
poderia, em princípio, ser muito similar àquela 
proveniente de um buraco negro, diferenças da 
dinâmica de fluxo na superfície estelar (H. Oli¬ 
vares et al. 2019, em preparação), fortes cam¬ 
pos magnéticos(Lobanov 2017), ou excesso de ra¬ 
diação no infraverme lho (Broderick k Narayan 
2006) permitiria distinguir um gravastar de um 
buraco negro. 

Em conjunto, os resultados aqui apresenta¬ 
dos fornecem uma nova maneira de estudar os 
espaços-tempos de objetos compactos e são com¬ 
plementares à detecção de ondas gravitacionais 
de coalescência de buracos negros de massa este¬ 
lar com LIGO/Virgo (Abbott et al. 2016). Nos¬ 
sas restrições aos desvios da geometria de Kerr 
dependem apenas da validade do princípio de 
equivalência e são agnósticas em relação à teo¬ 
ria da gravidade subjacente, mas podem ser usa¬ 
das para medir, com precisão cada vez melhor, os 
parâmetros da métrica de fundo. Por outro lado, 
as observações atuais de ondas gravitacionais da 
coalescência de estrelas investigam a dinâmica da 
teoria subjacente, mas não podem confiar na pos¬ 
sibilidade de medições múltiplas e repetidas da 
mesma fonte. 

Para sublinhar a complementaridade de ob¬ 
servações de ondas gravitacionais e eletro¬ 
magnéticas de buracos negros, notamos que uma 
característica básica dos buracos negros em RG é 
que seu tamanho escala linearmente com massa. 
Combinando nossas restrições no buraco negro 
supermassivo em M87 com aqueles nos buracos 
negros de massa estelar detectados via ondas gra¬ 
vitacionais, podemos inferir que essa propriedade 
permanece a mesma em mais de oito ordens de 
grandeza. Embora essa invariância seja uma pre¬ 
visão básica de RG, que é uma teoria livre de 
escala, ela não precisa ser satisfeita em outras te¬ 
orias que introduzem uma escala no campo gra¬ 
vitacional. 

Finalmente, o núcleo de rádio no M87 é bas¬ 
tante típico para poderosos jatos de rádio em ge¬ 
ral. Ela cai no plano fundamental da atividade 
dos buracos negros para núcleos de rádio (Falcke 
et al. 2004), conectando através de leis de es¬ 
cala simples as propriedades de rádio e raios X 
de candidatos a buracos negros de baixa lumi¬ 
nosidade em escalas de massa e taxa de acreção 
muito diferentes. Isso sugere que eles são alimen¬ 
tados por um objeto comum invariante à escala. 
Portanto, estabelecer a natureza do buraco negro 



para M87* também apóia o paradigma geral de 
que buracos negros são a fonte de energia para 
galáxias ativas. 

9 Conclusão e Perspectivas 


Montamos o EHT, um arranjo global VLBI 
operando em um comprimento de onda de 1,3 
mm e constituimos imagens de estruturas do ta¬ 
manho do horizonte de eventos em torno do can¬ 
didato a buracos negros supermassivo em M87. 
Usando vários métodos independentes de cali- 
bração, geração de imagens e análise, descobri¬ 
mos que a imagem é dominada por uma estru¬ 
tura em anel de 42 ± 3/ias de diâmetro, que é 
mais brilhante no sul. Esta estrutura tem uma 
depressão de brilho no centro com um contraste 
de > 10:1, que identificamos com a sombra do bu¬ 
raco negro. A comparação dos dados com uma 
extensa biblioteca de imagens sintéticas obtida 
a partir de simulações GRMHD cobrindo dife¬ 
rentes cenários físicos e condições do plasma re¬ 
vela que as características básicas de nossa ima¬ 
gem são relativamente independentes do modelo 
astrofísico detalhado. Isso nos permite obter 
uma estimativa para a massa do buraco negro 
de M = (6, 5 ± 0,7) x lO^M©. Com base em 
nossa modelagem e informações sobre o ângulo 
de inclinação, derivamos o sentido da rotação do 
buraco negro como sendo horário, ou seja, o ve¬ 
tor de giro do buraco negro aponta para longe 
de nós. O excesso de brilho na parte sul do anel 
de emissão é explicado como uma irradiação re- 
lativística do material girando no sentido horário 
como visto pelo observador, isto é, a parte in¬ 
ferior da região de emissão está se movendo em 
direção ao observador. 

Observações futuras e análises adicionais tes¬ 
tarão a estabilidade, a forma e a profundidade da 
sombra com mais precisão. Uma de suas princi¬ 
pais características é que deve permanecer am¬ 
plamente constante com o tempo, uma vez que 
não se espera que a massa de M87* mude de ma¬ 
neira mensurável em escalas de tempo humanas. 
A análise polarimétrica das imagens, que repor¬ 
taremos no futuro, fornecerá informações sobre a 
taxa de acreção via rotação de Earaday (Bower 
et al. 2003; Marrone et al. 2007; Kuo et al. 2014; 
Moscibrodzka et al. 2017) e no fluxo magnético. 
Imagens de resolução mais alta podem ser obtidas 
com um comprimento de onda menor, ou seja, 0,8 
mm (345 GHz), adicionando mais telescópios e, 
em um futuro mais distante, com interferometria 
dependente do espaço (Kardashev et al. 2014; 
Eish et al. 2019; Palumbo et al. 2019; E. Roelofs 
et al. 2019b, em preparação). 


Outra fonte primária do EHT, Sgr*, tem uma 
massa precisamente medida com três ordens de 
magnitude menor que a de M87*, com escalas 
de tempo dinâmicas de minutos em vez de dias. 
Observar a sombra da Sgr A* exigirá a contabi¬ 
lização dessa variabilidade e mitigação dos efeitos 
de espalhamento causados pelo meio interestelar 
(Johnson 2016; Lu et al. 2016; Bouman et al. 
2018). Análises de dependentes de tempo mas 
não através de imagens podem ser usadas para 
rastrear potencialmente o movimento de emissão 
de matéria perto do buraco negro, conforme rela¬ 
tado recentemente através de observações inter- 
ferométricas no infravermelho próximo (Gravity 
Collaboration et al. 2018b). Tais observações for¬ 
necem testes independentes e sondagens de RG 
em ainda outra escala de massa (Broderick & 
Loeb 2005; Doeleman et al. 2009b; Roelofs et 
al. 2017; Medeiros et al. 2017). 

Em conclusão, mostramos que os estudos dire¬ 
tos da sombra do horizonte de eventos de candi¬ 
datos a buracos negros supermassivos são agora 
possíveis através de ondas eletromagnéticas, 
então transformando este limite difícil de vi¬ 
sualizar, originário de um conceito matemático 
em uma entidade física que pode ser estudada 
e testada através de repetidas observações as¬ 
tronômicas. 
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